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Q^ ■ Abstract 



Study of the MSSM in large tan/3 regime has to include correlations between the 
_^ , constraints presented by the low energy values of the b quark mass and BR(6 — > 57). 

QJ , Both quantities receive SUSY contributions enhanced by tan/3 and have a major impact 

on the MSSM analysis. Here we summarize the results of such a study constrained 
by a complete SO (10) model. The dominant effects to the analysis come from third 
p\ ' generation Yukawa couplings and a GUT threshold to a^. We show that a small 

C^ ■ negative GUT correction to ag accommodates as{Mz) < 0.118 and 6m,^^^^ > 0. The 

latter quantity being positive then opens up two options to fit the measured rate for 
6 — > 57. The two distinct fits differ by the overall sign of the amplitude for this process. 
They work equally well in complementary regions of the allowed SUSY parameter space. 
We show plots of the partial contributions to the coefficient Cj{Mz) in the {mo, M1/2) 
plane in each of these best fits. We conclude that an attractive SO(10)-derived regime 
of the MSSM remains a viable option. 



PACS numbers: 12.15.Ff, 12.15.Hh, 12.60.Jv 

*0n leave of absence from the Dept. of Theoretical Physics, Comenius Univ., Bratislava, Slovakia; 
blazek@gluon2.physics.indiana.edu 



raby@pacific.nips.ohio-state.edu 



1 Introduction 

The inclusive b ^ s'j decay has attracted a lot of attention in the past and will undoubtfully 
attract at least as much attention in the future. On one side, it is an observed flavor changing 
neutral current (FCNC) process with the measured rate BR{b -^ S7) = (2.32±0.57±0.35) x 
10"^ announced by CLEO0], and BR{b -^ 57) = (3.38 ± 0.74 ± 0.85) x 10"^ found most 
recently by the ALEPH Coll. , and we can expect that the statistical uncertainty of these 
measurements will be reduced in the near future. On the other side, the Standard Model 
(SM) prediction with the next-to-leading order QCD corrections included has been calculated 
as BR{b — > s'-f) = (3.48 ± 0.31) x 10~^ §]. However, since forbidden at tree level it lets 
the SM compete with the loop contributions from any new physics. Theoretical analyses of 
this process, which assume minimal supersymmetric (SUSY) extension of the SM and are 
constrained by unification, naturally fall into two categories. Either tan/5 (the ratio of the 
Higgs vacuum expectation values <H^> / <H^> ) is considered to be low {i.e., close to 1) 
or this parameter is large, of the order 50. 

In the first case, both < i7° > and < H^ > are of the order of the electroweak scale and 
the hierarchy Mf,, M^ <C Mt is triggered by the smallness of the b and r Yukawa couplings Xf, 
and At-. In effect, this regime of the Minimal Supersymmetric Standard Model (MSSM) has 
several simple features: radiative electroweak symmetry breaking is rather straightforward to 
impose, and the terms suppressed by Af, and A,- can be neglected both in the renormalization 
group equations (RGEs) and in the analysis at the Z scale. A drawback, if we can say so, 
is that when considering the MSSM as a low energy effective theory, there is no appealing 
unification dynamics requiring the hierarchy A;,, A,- <C Atop- 

The large tan/3 case, on the other hand, is very attractive when considered from the 
unification point of view. Using the MSSM RGEs one finds that this case offers an amazingly 
simple option A;, = A,- = Atop at the unification scale, consistent with the minimal 5*0(10) 
grand unified theories (GUTs). The third generation mass hierarchy is, however, explained 
by < H^ > ~ 3-4GeV, the value which is much less than the electroweak scale - the generic 
scale of the other dimensionful parameters, including < H^ >. That makes the MSSM 
analysis with large tan/5 challenging in more than just the symmetry breaking sector. For 
example, there are fermion masses which are proportional to <H^> at tree level, and there 
is no symmetry which would guarantee the same suppression in their radiative corrections. 
Contributions at one- loop may be and indeed are proportional to < H^ >. The same type 
of corrections are then also induced to some CKM matrix elements, e.g. Vcb or Vub- As a 
result, the tan/? enhanced corrections to such observables must be included into the MSSM 
analysis. Inclusive 6 ^ 57 decay represents an even more sensitive probe if tan/? is large. As 
an FCNC process emerging only at a loop level its SM contribution turns out to be suppressed 
by < if ° >. On the other hand, the chargino-squark contributions are proportional to < i/° > 
and may dominate the whole process. Thus the MSSM analysis becomes more powerful (and 
restrictive) than in the case with low tan/5. 

In this work, we present the results of such a constrained global analysis of the MSSM 
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with large tan/3 and pay special attention to the structure of the b ^ s'y partial amplitudes. 
Our actual investigation has been performed within a complete 5*0 (10) model proposed by 
Lucas and RabyQ. The model (called model 4c), was analyzed in detail in 0] and a very 
good agreement with low energy data (including fermion masses and mixings) was found in 
a large region of the SUSY parameter space. Although the primary goal of the analysis was 
to discriminate among different models of fermion masses, this model turned out to work so 
well, that the main constraints were the b quark mass and the BR{b —>■ sj). As a result we 
can regard the work as testing the MSSM with large tan/3, with the inclusion of a particular 
set of 3x3 non-diagonal Yukawa matrices at the GUT scale. Thus we think that this letter 
naturally fits into the mosaic of many previous works on 6 — i> 57 in the MSSM. These works 
studied the 6 ^ 57 process either with low (or moderate^, i.e. less than about 10) tan/5 
onlyP] or, if the large tan/? regime was also considered 0, then some of the important 
ingredients of our procedure were not taken into account. Most notably, we improve the 
previous studies (i) by taking into account the tan/3 enhanced SUSY threshold corrections 
to fermion masses and mixings, and (ii) by introducing global analysis instead of scatter 
plots of random points in SUSY parameter space. Due to the latter improvement we can 
present the contour plots of the interesting quantities from the best fits and observe how 
these quantities are correlated. Additional motivation for this paper came from the study 
P|. This work has also introduced global analysis, with one-loop SUSY threshold corrections 
properly included, and studied both the low and large tan/3 regimes. However, this analysis 
assumed strict gauge unification with no threshold correction to a^ at the GUT scale. As a 
result, it did not concentrate on the regime with rather low values of as{Mz) < 0.118 and 
positive SUSY corrections to mb. Therefore, that work does not test the subspace in SUSY 
parameter space which we study in this analysis^]. 

Thus our main result is that in the regime of SUSY parameter space where tan/3 is large 
and the b quark mass receives positive corrections we find the analysis consistent with t — b—r 
unification and the observed BR(6 -^ S7), in contradiction to previous studies. Dominant 
effects come from the third generation Yukawa couplings and the introduction of 63, a GUT 
threshold to a^. In the best fits, small negative €3 is correlated with lower ^^(M^), and that 
in effect decreases irib. We also find that the contribution of the inter-generational ci-ti 
squark mixing to the b ^ s'j amplitude (tan/3 enhanced compared to the SM contribution) 
is numerically significant. It is, however, model dependent. Note that it is neglected in 
the approximation of Barbieri-Giudice which has been used in many of the succeeding 
studies. 



^with no reference to the unification of the Yukawa couplings 

^We have, in fact, observed identical features in the same subspace of the SUSY parameter space and 
have come to the same conclusions as in [^ that the subspace where 5mf^^^ < appears to be excluded 
— see the discussion on negative /i parameter in section 3 and in B. (In our conventions, /i < implies 
negative SUSY corrections to the b quark mass and purely constructive interference among all leading MSSM 
contributions to the decay b — > s"f.) 



The main steps of our procedure are described in section 2. For completeness, the per- 
formance of model 4c in the global analysis is presented in section 3. Section 4 contains the 
results of this work. It starts with comment on model sensitivity and includes contour plots 
in the (mo, M1/2) plane of constant BR{b — » 57) and various contributions to this process, 
extracted from our best fits. We show how the destructive interference among the b —* sj 
partial amplitudes provides for two distinct fits, each of a very low x^ iii rather complemen- 
tary regions of the parameter space. We also discuss phenomenological implications of our 
findings. Finally, conclusions in section 5 contain a brief summary of this work. 



2 Global Analysis 

Details of our numerical analysis are described in 0. Here we summarize the main steps 
of our procedure relevant for this letter. We perform a global analysis in a strict top-down 
approach. We start with the following initial parameters: the scale of new physics Mg, 
unified gauge coupling aciMc) 0, one-loop GUT threshold correction €3 to 0^3 (Mg), and 
A: the 33 element common to all three Yukawa matrices at Mq- We assume supergravity 
induced SUSY breaking and neglect the effects of running between the Planck scale and the 
GUT scale. The parameters of the SUSY sector, which are introduced at Mq, include a 
common gaugino mass M1/2, common scalar mass niQ of squark and slepton mass matrices, 
scalar Higgs mass parameters friH^^ and rnn^ and a universal dimensionful trilinear coupling 
Aq. For simplicity, the n parameter and its SUSY breaking bilinear partner B are introduced 
at the Z scale, since they are renormalized multiplicatively and do not enter the RGEs of 
the other parameters. In the actual model 4c analysis, the structure of the Yukawa matrices 
at Mg is built up using six more dimensionless parameters. As explained in section 4, these 
are small numbers and their effects decouple from the observables related to the gauge and 
SUSY sectors and from the masses of the third generation fermions. 

We use the two-loop RGEs for the gauge and Yukawa couplings, and the one-loop RGEs 
for the MSSM dimensionful parameters to run down to the Z scale. At selected points, 
we check that the full two-loop RGEs of the MSSM ||10| yield the same results. At the Z 
scale, the MSSM is matched to the effective theory consisting of QCD and electromagnetism, 
leaving out the SM as an effective theory. Within the MSSM, we implement the effective 
potential method of [|IT] (at one loop, with all one-loop threshold effects included) to obtain 
the fit values of v and tan/3. These quantities are not directly restricted. (Neither is tan/3 
among the free initial parameters.) Electroweak symmetry breaking is established implicitly 
in the process of the x^ minimization by the evaluation of the fermion and gauge boson 
masses. The latter are calculated with the full MSSM one-loop corrections included. We 
also compute all the threshold corrections, proportional to tan/3, to fermion masses and 
CKM matrix elements, following [Q. Masses of the Higgs particles are calculated at the 



^Mq is defined as the scale where the gauge couphngs ai and a2 are exactly equal within the one-loop 
GUT threshold corrections. By ac we actually mean the value ai{MG) = ct2{MG) 



Observable 


Central 


a 


Observable 


Central 


a 






value 








value 




1. 


Mz 


91.186 


0.46 


11. 


Mb-Mc 


3.4 


0.2 


2. 


Mw 


80.356 


0.40 


12. 


nis 


180 


50 


3. 


G, 


1.166-10-^ 


1.2- 10-^ 


13. 


rud/ms 


0.05 


0.015 


4. 


a-' 


137.04 


0.69 


14. 


Q~' 


0.00203 


0.00020 


5. 


as{Mz) 


0.118 


0.005 


15. 
16. 


Me 


105.66 
0.5110 


0.53 
0.0026 


6. 

7. 


Mt 
mb{Mb) 


175.0 
4.26 


6.0 
0.11 


17. 


Vus 


0.2205 


0.0026 


8. 


Mr 


1.777 


0.0089 


18. 
19. 


Vub/V,b 


0.0392 
0.08 


0.003 
0.02 


9. 


Pnew 


-0.6 • 10-3 


2.6 ■ 10-3 


10. 


B{h -^ S7) 


2.32-10-^ 


0.92 ■ 10"^ 


20. 


Bk 


0.8 


0.1 



Table 1: Experimental observables of the global analysis. 



two- loop level as in [|r^ , and masses of the SUSY particles are left at their respective tree 
level values. The tree level squark and slepton masses are constrained to be greater than 
30 GeV. This is below the experimental limit but when they are that light, we count on 
substantial enhancements at one loop |jl3[. Chargino and neutralino masses receive very 
small one-loop corrections and so we restrict their tree masses by the respective LEP limits. 
As we match the MSSM to the effective theory below Mz-, the threshold corrections to the 
gauge couplings are calculated following |T^, with the exception that the SUSY vertex and 
box corrections to Ar are neglected. In our approach, Ar serves to derive the theoretical 
value of G^. When below the Z scale, we use the three-loop QCD and one-loop QED RGEs 
Tsj to evaluate the quark and lepton masses. 

Our y^ function is calculated based on the low energy data (observables and their corre- 
sponding errors) listed in table |l]. Our analysis was originally designed to test GUT models, 
so from the point of view of testing the MSSM the data in table are divided into two 
groups, corresponding to observables 1-10, and 11-20 respectively. Five observables in the 
gauge sector (M^, My/, G^, a and a^), masses of the third generation fermions, pnew Q and 
BR{h —>■ s'-f) are typically chosen to test the MSSM constrained by unification. The other 
ten observables corresponding to six light fermion masses and four independent parameters 
of the CKM matrix have been included in the analysis but they do not significantly affect 
our results (see introduction to section 4 where we discuss model dependence). 

Note that seven out of the twenty observables have the estimated theoretical uncertainties 
dominating over the experimental ones (their respective ex's are underlined in table |l|). These 
theoretical uncertainties represent conservative estimates of the errors (0.5% for six out of 
the seven, and 1% for G^ to compensate for SUSY boxes and vertices which are neglected 
in the computation of Ar) generated by our numerical procedure. In addition, note that we 



*This is the contribution of physics beyond the SM to the p parameter. 



have also introduced a conservative error on as{Mz) [|T^ and added the CLEO errors for the 
BR{b — i> 57) hnearly. 

2.1 Calculation of 6 ^ 57 



The MSSM amphtude for the transition 6 — i> 57 is calculated following [T^ at the threshold 
Mz- The effective Hamiltonian method, summarized in [|18|, is used below Mz- In particular, 
the amplitude at Mz is matched to the Wilson coefficient Cj{Mz) in the effective hamiltonian 



H, 



eff 



AC ^ 



Following the conventions of ref . ||T3 the magnetic dipole operator reads 



O, = ^^m,,{si^a^'hn)F^, , (2) 



with a^^ = i/2 [7/^,7;/]. The branching ratio is computed from the formula 



I T /^^ T /" I 2 /"' 

BRih ^ S7) = ^-^^ Tir^TTTT 1^7^^ (/^fc) I' BR(h -. cev) , (3) 



where BR(h — > cez/) = 0.104 [19|, a = 1/132.5 and the phase-space function for the semilep- 



tonic decay g{z) = 1 — 9>z + 82 — z — 2Az \ogz. The effective coefficient 

Cf^(^) = rj'iCriMz) + ^ (r^M - n'i) C^{Mz) + C^iMz) j2h^v''' (4) 

comes from operator mixing in the leading log approximation, with rj = as{Mz)/<ys{fJ')- The 
numbers hi and a^ are given in |]T3. For the CKM matrix elements and quark masses, the 



fit values of the model 4c analysis are consistently used in the formulas above. The best fits 
yield their values very close to the ones quoted by Particle Data Group (PDG) [|1^. The 
values of as{Mz) are also taken from the GUT analysis. (The best fit values which are used 



to calculate rj in eq.(|D can be read out from the figures discussed in section ^^ .) 

In our analysis, we fix the low energy scale to be fib = 4.7GeV and do not study the 
scale dependence of the result. The scale dependence will be reduced once the complete 
next-to-leading order calculation within the MSSM will be known. Because of the significant 
SUSY contributions to this process in large tan/3 regime, we use the leading order calculation 
in our fits and only comment (in section [4.5|) on possible changes which may result from the 
next-to-leading order calculation in the full MSSM. 

^ Since we do not choose some fixed values, we observe a moderate dependence of the BR{b — > 57), within 
10%, on the particular set of the computed theoretical values for these quantities. For the most part, it is 
due to the function g(z) which changes rather fast in the vicinity of zwO.S. 



3 Model 4c Best Fits 

The analysis, as described above, was used to test simple SO (10) models [^. It was found 
0, ^ that one of the models, called model 4c 0], yields very good fits in a large portion 
of the allowed SUSY parameter space. The quality of these fits is presented in figures la-c. 
The figures show the contour plots of the minimum x^ in the {rriQ, M1/2) plane, for three 
different fixed values of the parameter fi{Mz) = 80, 160, 240 GeV. All initial parameters 
other than {mo, M1/2, ft} were subject to minimization. 

Note that as /i increases, the quality of the fit gets worse. This is understood from the 
form of the SUSY corrections to fermion masses and mixings. These corrections increase with 
fi and as fi gets larger they can only be kept under control by larger squark masses. For this 
reason the contour lines of constant x^ move towards larger values of mo for M1/2 < 400GeV 
in figures lb and Ic. Varying fj, freely actually results in its approaching the lowest possible 
value. This lower bound on ^{Mz) is determined by the chargino mass limit from direct 
searches and is correlated with Mi/2- When the value of fi is fixed, as in figures la-c, the 
chargino mass limit then sets a sharp lower bound on M1/2, which is explicitly visible in each 
of the figures. Plots in figures la-c were constructed under the assumption that the chargino 
mass limit was 65 GeV. 

Because the chargino mass limit has been raised at LEP2, and because the optimization 
ends up with low values of /i, in the rest of the analysis presented in this paper fi{Mz) has 
been fixed to llOGeV. At this value of /x, the lightest chargino mass turns out to be about 
lOOGeV for M1/2 > 340GeV, and slowly drops down to about 85GeV for M1/2 = 200GeV. [] 
Figures 2a and 2b show explicitly that the structure observed in figures la-c originates from 
the two distinct fits corresponding to two separate minima of the global analysis. The fits 
are primarily distinguished by the sign of the effective 6 — *> S7 decay amplitude, or in other 
words, by the sign of the C7 coefficient of the effective Hamiltonian for the low energy FCNC 
processes. The two options are available because of the destructive chargino interference 
among the partial amplitudes to this process. A more detailed discussion of these results 
is actually the subject of section 4. As can be seen from the contour lines corresponding 
to x^ =0.4 and 0.3, the minimum in the {mo,Mi/2) plane is quite shallow. For this reason 
and because the optimization converges very slowly we do not specify the exact point with 
minimum x^ in the figures. We are also aware that such low x^ values most likely result from 
an overestimate of the theoretical uncertainties. That suggests that our best fit x^'s should 
be used more in the sense of figures of merit than in a rigorous statistical fit evaluation. 

We do not show results for negative values of /i. In this case, the SUSY corrections to m?, 



^Tlic figures are taken from m. These are slightly modified compared to the figm'es in pTl due to a sign 
error found later in the numerical code. 

^The current limit is m^- > 87-90GeV (dependent on tan/3 and the rest of the SUSY spectrum) from 

the LEP2 run at ^/s = 183GeV. [^ We would need to increase our value of /i by a few GeV in order to get 
over this limit in the region where M1/2 = 200-225GeV. Such a change would, however, be insignificant for 
the rest of our results. 



are negative (which is rather a welcome feature in connection with strict b — t unification). 
However, the chargino contribution to b ^ s'j interferes constructively with the already 
large enough SM and charged Higgs contributions. As a result, the fits get much worse, with 
X^ well above 10 per 3 d.o.f., and that disfavors this region of the SUSY parameter space. 
Similar observations were also made in PI. 



4 MSSM Analysis 

4.1 Model Dependence 

Recall that a general, model independent analysis of the MSSM constrained by unifica- 
tion traditionally assumes exact gauge and Yukawa coupling unification and some degree 
of universality among the SUSY mass parameters. The lighter generation fermion masses 
and the CKM matrix elements are taken over from experiment, and in the SUSY sector, 
only the left-right mixings of the stops, sbottoms and staus are usually considered, with the 
inter-generational mixings left out. 

The most significant model dependent feature of the analysis presented in this letter 
is the introduction of 63, the GUT threshold to Og- We have introduced 63 as an initial 
parameter which is free to vary within ±6%. The contour plots of constant €3 resulting from 
the best fits are shown in the figures discussed in section ^^. As one can see, €3 tends to 
be negative. Lucas and Raby showed that such negative threshold corrections are consistent 
with the complete SO{10) formulation of model 4c 0|. Negative €3 allows as{Mz) to go 
below 0.120. That in turn is a welcome feature for the b — t unification if one studies the 
SUSY parameter subspace in which mi, receives positive SUSY corrections. There is more 
discussion on these effects in section |4.6|. €3 is the only GUT threshold introduced in this 



study. 

For the Yukawa matrices, the exact equality of the 33 elements is assumed. The remaining 
Yukawa entries are, of course, dependent on specific properties of model 4c. However, they 
are small and decouple from the MSSM RGEs for the gauge and third generation Yukawa 
couplings, as well as for the diagonal SUSY mass parameters. Thus they have no effect on the 
calculation of the Z-scale values for the first nine observables in table |l|. Only the branching 
ratio BR{b -^ 57) is affected by some of these entries. This dependence comes dominantly 
from the diagram in figure ^. The partial contribution to b -^ s'-f is proportional to the 
fiavor changing CL-t^ squark mixing in this case |2^ and is tan/3 enhanced which makes it 



non-negligable. The mixing is completely induced by the off-diagonal entries of the Yukawa 
matrices in the RG evolution, since we assume universal squark masses at the GUT scale. In 
addition, there is no significant pull in the model 4c best x^'s from the light fermion masses 
and CKM mixing elements [^. Instead, the dominant pulls are imposed by as{Mz), mf,{Mi,) 
and other observables which also enter the MSSM analysis. Typically, about 60-80% of the 
total x^ value comes from the first ten observables of table |l]. That means that it is these 
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traditional MSSM observables which drive the model 4c optimization procedure when it gets 
close to its minima and that at the same time the Yukawa sector of model 4c works indeed 
very well and does not bias the optimization significantly. We conclude that our results 
presented in this section are not sensitive to the structure of the Yukawa matrices except for 
the model dependent 23 mixing which is significant for the BR(6 —>■ s'j). 

4.2 BR{b -^ 57) in MSSM with large tan/3 

The SM, charged Higgs, and chargino diagrams contributing to the amplitude for b —>■ s'-f are 
shown in figures ^-e. These are the dominant contributions in the best fits of our analysis. 
Note that the three chargino contributions in figures ^-e are enhanced by tan/3. In our code, 
we also take into account the gluino and neutralino diagrams as well as the full chargino 
contribution including the pieces not enhanced by tan/3, although their contributions are 
numerically insignificant. 

To understand our results which will be presented below, we would like to look first at 
an estimate what to expect from the SUSY contribution to C^lMz). As a starting point we 
acknowledge the fact that the SM contribution gives a rough agreement with the measured 
rate and define the ratios for Higgs and chargino contributions 

rW = cfVc^f*'\ (5) 

r(^) = cP/dr'''\ (6) 

where all quantities are evaluated at Mz- Equation (^) then reads 

where we used C2{Mz) = 1 and, for simplicity, neglected the mixing with the chromomag- 
netic operator, proportional to Cg, which turns out to be by about a factor 5-10 less than 
the other two terms. 

Since Cj ^ (approximately equal to 77 23 C*} +J2 hii]^' ) would yield about the right 
value of the BR{b -^ sj) we deduce that either 

for C\ ~ —C\ >^J\ Por the last estimate, the numerical results C\ = —0.190, 

7723 = 0.679, and I]i=i hii]""- = —0.168 were used — computed for as{Mz) = 0.118. 

The charged Higgs contribution always interferes constructively with the SM contribution 
0. Typically, we get 

0<^W<1.3^ (10) 



depending on the mass of the H~. In the first case, especially if r^^^ and r^^^ are non- 
negligible, eq.(|D means that the chargino part must interfere destructivelyj^with the SM and 
charged Higgs contributions, practically cancelling the latter. The enhancement by tan/5 of 
the chargino contribution has to be compensated for by rather large masses of the sparticles 
in the diagrams in figures 0c-e. In the second case, described by eq.(^, large destructive 
chargino interference is required to outweigh the combined SM and H~ contributions and to 
flip the overall sign of the amplitude. Quite amazingly, it is not so difficult to arrange (see 
also [^) since the chargino contribution is the only one enhanced by large tan/?. However, 



large sparticle masses obviously suppress the effect. The lesson is that we can expect the 
two cases to work in a complementary SUSY parameter space and have sign Cj =(-) 

sign C7 for the best fits in the region with large (low) uiq and/or M1/2, respectively. 

These expectations are indeed realized in the best fits of model 4c. Figure 2a shows 
the best fits in the case when the chargino contribution roughly cancels the charged Higgs 
contribution, while the best fits of figure 2b correspond to the case when the chargino piece 
truly dominates and reverses the sign of the overall amplitude. As anticipated, these two 
cases work in complementary regions of parameter space. 

Figures ^a and ^d show how well the fits describe the measured value of the BR{b —>■ s'j). 
In these figures (and similarly in the following ones) we show the contour lines of constant x^ 
in the background for reference. As can be seen, the BR{b — > S7) indeed presents a major 
constraint since whenever the x^ values go up, the agreement with the observed b ^ sj 
decay rate gets worse. Figures |]a-b and ^a-b show the contour plots of constant r^^^ and 
f-iC) in each of these two cases. One can compare the numerical results in these figures with 
the approximate relations (||) and (H). Note also the validity of relation ([lOD in each case. 
Finally note that due to large tan/5 the effects of SUSY decoupling start showing up only 
for mo > 2TeV, i. e. outside the SUSY space studied in these figures. 

4.3 Discussion of the results for BR{b -^ 57) and phenomenological 

implications 

There is one striking feature which is common to both cases. It is that both fits would like to 
have the BR{h — > 57) below rather than above the current experimental value 2.32 x 10"'*. 
In the first case, the tan/3 enhanced chargino contribution tends to be too large when 
going against the charged Higgs contribution, since the latter is not tan/? enhanced. The 
fit clearly favors as large Higgs contribution as possible with r^^'^ reaching its maximum 
(see fig.|^). A phenomenological consequence of this observation is that the charged Higgs 
(and then also the whole Higgs sector) tends to be as light as possible in this case. We get, 

^Since eq's (M) and (ph are valid only approximately, the case 1 + r*^^) + A'-^^ « 1 in principle also allows 
for a constructive interference < r^'-^' , r*-^-* ^ 1 in the region in parameter space where rriH- and sparticle 
masses are large. This option, however, does not result from our best fits, as already mentioned in the 
discussion on negative ^ parameter in the previous section. 



for instance, the best fit value of the pseudoscalar mass niA < lOOGeV everywhere in the 
(mo, M1/2) plane. 

In the second case, when the chargino part overshoots the combined SM and H~ con- 
tributions we observe different effects in the regions with M1/2 below and above (roughly) 
SOOGeV. For larger values of M1/2 the chargino contribution clearly tends to be not big 
enough. As a result we might expect to see only very low values of r^^^ in this region — 
complementary to the large values in fig.^a — and a very heavy Higgs sector. However, the 
best fit value of r^^^ varies quite a bit indicating that the charged Higgs mass does not stay 
at some very large value. That is related to the observation [^ that one cannot have good 
fits with the Higgs sector much heavier than squarks. When M1/2 gets below SOOGeV the 
b ^ s'j decay rate is no longer a strong constraint and two separate minima can be found in 
the course of the optimization. The two fits work equally well: x^ i^i each case stays below 1 
per 3dof. The minima differ by the best fit values of the Higgs masses: one minimum corre- 
sponds to rriho and rriA at the experimental lower limit (set to 65GeV in our analysis) while 
these masses gradually rise in the second "valley" up to 700GeV. When crossing the region 
with Ml/2 ~ SOOGeV towards larger M1/2, the first "valley" vanishes and the optimization 
slides down to the second minimum because of the BR{b —>■ sj) . In the allowed corner with 
rrio < 700GeV, the two "valleys" approach each other and finally coincide. The effect of the 
doubled minima is indicated in figures §3 and Pd with the solid black (dashed gray) contour 
lines corresponding to the heavier (lighter) Higgs sector^. 

In summary, if the future experimental analysis confirms the discrepancy between the 
CLEO measured value and the NLO SM calculation 0, the MSSM with large tan/5 could be 
the solution. Similarly, if the NLO calculation is completed for the MSSM, and if it turns out 
to increase the LO result as occurred for the SM, then the large tan/3 regime will apparently 
have no problem fitting the 6 — > S7 rate exactly. That seems to be in contrast with the fits 
in the low tan/? regime of ref.0, which get below the SM value only for M1/2 < 200GeV and 
small tuq. 

4.4 Role of c-t mixing in our results for BR{h -^ 57) 

It is interesting to note the significance of the inter-generational squark mixing in fig.^. In 
analogy to eq's (|^), (|^) we define 

^(C23) ^ C'P^^) / ^(5M) ^ (11) 

where C) is the ci-ti mixing contribution to the coefficient C-j at the scale Mz- We show 
the contour plots of the constant r^*^^^) in figures |^ and |^. While the dominant chargino 
contribution is that proportional to the t^-tij mixing (fig.Sc), C\ becomes important 
because of the destructive character of the interference among the partial amplitudes. As one 
can see, the ci-ti mixing term can be comparable with the SM and H~ contributions. More 



'The same holds in figures ^b and ^. 
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importantly, it always interferes constructively with them. In the case when sign C} = 

+ sign C7 , this term helps to counterbalance the large contribution of the left-right stop 
mixing. In the complementary case, when the chargino contribution overturns the sign of 
the net amplitude the CL-ti mixing makes this flip more difficult to happen. As a result, 
it has different consequences for the two fits in figures 2a and 2b, especially important in 
the region (mg < lOOOGeV, M1/2 ~ 350GeV) where the fits start getting worse. In figure 
2a, it improves the fit for lower values of (mo, M1/2) where the contribution from the stop 
mixing alone would otherwise overwhelm the sum of the SM and H~ diagrams. On the 
contrary, it worsens the fit in figure 2b in the same parameter subspace. These observations 
are, however, model dependent as has been explained in the introduction to this section, 
and their validity relies on the boundary conditions assumed at the GUT scale. A general 
study of the limits imposed by the inter-generational mixings can be found in the review ^| 



and the references therein. To make a connection with the general approach we note that 
typically we get 



(m 



2 \23 



iSu)LL ^ ., 2X2/ 2^33 11/2 ~ -(0-01-0.02) (12) 

[(m^)ii(m^)fJV2 

at the Z scale, where {'mfj)^ is the 3x3 rriQ squark mass matrix after the unitary rotations 
which diagonalize the fermionic sector are performed in the squark sector too ( — sandwich- 
ing by Vff's in this particular case). This value is to be compared with the fiavor-changing 
effects originating in the CKM matrix (figures |^a-d) with \Vts\ ~ \Vcb\ ~ 0.04. 

4.5 Effects of NLO QCD Corrections to BR(6 -^ 57) 

Finally, we would like to comment on the next-to-leading order (NLO) QCD corrections. 



These have been completed only for the SM 0]. In ref. [^ it was shown that the NLO 
matching at the high scale gives a non-negligible contribution to C7. In particular, for the 
SM we get 

= -0.312 + 0.017(-0.883-0.111 + ---) , (13) 

where Cj (jJb) is given in terms of the LO matching coefficients (our eq.(^) and for the 
evaluation in the second line as{Mz) = 0.118 was assumed. Hence within the SM the NLO 
matching corrections C7 (Mz) and Cg (Mz) alone would change the rate by about 10%. 
The final NLO corrected result increases the LO value obtained at fib — Mb by about 20%. 
(In other words, it effectively lowers the scale fib of the LO calculations to about 0.6 Mb as 
obtained also in ref.|p.) In the MSSM, the NLO matching remains to be calculated. Clearly, 
with large tan/? the chargino diagrams will likely dominate and their contribution may cause 
a different effect than in the SM. For that reason we have left our results without applying 
the NLO corrections and will just sketch at this point what effects the higher order matching 
corrections may have in our best fits. 
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The case when sign Cj = + sign C-j- is much hke the SM case. The effect of the 

unknown C^ \Mz) and Cg (Mz) can be readily estimated from the complete 

form of eq.([T3|) quoted in [Q. If the NLO matching followed the cancellation among the 
charged Higgs and chargino contributions observed in the LO (eq.®), our results for the 
BR{h -^ 57) would also be about 20% larger, similar to the SM calculation. We can see from 
figure ^a that this would increase the parameter space in this case. Effectively, the NLO terms 
would add to the LO contributions of W~ and H~ , enabling the tan/3 enhanced chargino 
contribution to be greater in magnitude than what is allowed by eq. (H) — a welcome feature 
for lighter SUSY spectra. Clearly, if the NLO matching terms remain negative and get even 
larger in magnitude than in the SM case, the fit will improve substantially in the region with 
low ttlq and M1/2 (and may get somehow worse in the upper left corner of the (mo, iVfi/2) 
plane). On the other hand, we have checked that our LO results remain unaltered for 
C^ H j ^ —Q^^ , A; = 7, 8: for such matching terms all NLO corrections practically 

sum up to zero. Only if the NLO matching terms turn out to be very large and positive, the 
fit will be forced to retreat significantly towards larger M1/2 values. 

The second case, sign C7 = — sign C7 , is obviously more sensitive to the un- 

known NLO matching. Despite the sensitivity our sample calculations indicate that a drop 
in the value of the BR{h -^ 57) is fairly common and more likely to happen than an enhance- 
ment. For instance, we get BR{h — > S7) x 10'^ = 0.98 (2.31) in the complete NLO calculation 
assuming the NLO matching terms are Cf, = +1 (—7) x C^ , A; = 7, 8, and 

keeping cj"^ and Cf^ at the values which yield BR{h -^ 37) x 10^ = 2.31 at the LO. The 
reduction in the rate is correlated with the opposite signs of C-j and C% as compared to the 
SM case. However, it still holds that the NLO correction would be effectively taken into 
account by lowering the scale in the LO calculation. For the region with M1/2 > 300GeV, 
where the fit gradually gets worse in this case, it means that the 6 — *> 57 decay is likely 
more constraining than it appears in the analysis based on the LO calculation. As indicated 
above, the NLO corrections will improve the fit only if the NLO matching terms are at least 
seven times larger and opposite in sign compared to the SM NLO matching terms. 

4.6 (5mf^^^, as{Mz) and 63 in the Best Fits 

For the presented results, it has been important to have a destructive interference among 
the partial contributions to 6 — ;> 57. With the universal boundary conditions at the GUT 
scale however, that can be arranged only for a specific sign of the /x parameter: /i > in our 
conventions. That in turn correlates with the positive sign of the SUSY corrections to the h 



quark mass ||26[ . This fact implies strong constraints on the SUSY parameter space from tti^ 
since 5mf^^^ gets a dominant contribution from the tan/3 enhanced gluino exchange, which 
tends to be very large ||^ and offsets the agreement between the low energy value of irti, and 
the h-T unification at the GUT scale. The gluino correction can be explicitly suppressed by 
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heavy squarks and by the chargino correction which enters with the opposite sign^ It can 
also be reduced by lower values of as{Mz) because of the strong couplings in the vertices of 
the gluino - 6-squark diagram. The contour plots of the total SUSY correction to mh{Mz) in 
the best fits are shown in figures |a-b where one can see that these effects are quite effective 
in reducing the value of Smb. 

Note that the impact of large positive 5mf^^^ is reduced by a lower value of as{Mz) 
also indirectly — due to the RG evolution of the current h mass from Mz down to M^. The 
effect is enhanced if the current mass nib is converted to the perturbative pole mass Mf, in a 
top-down analysis. For example, the difference between as{Mz) being 0.115 and 0.121 leads 
to about 5% difference in M^, if the same value of mi,{Mz) is assumed in each case. This is a 
significant effect. In our analysis, which assumes larger uncertainty for a^ than for mf,{Mf,), 
it pushes as{Mz) down. The effect can be seen in figures ^-b. Note that one can trade 
lower values of as for higher values of rrib provided a smaller uncertainty o"(as) is assumed 

Rather low values of as{Mz) are, however, difficult to obtain from the exact gauge cou- 
pling unification, traditionally assumed in an MSSM analysis constrained by unification. 
In our analysis, we assumed that there is a few per cent correction to the gauge coupling 
unification generated by the spread in masses of heavy states integrated out at the GUT 
scale. It turned out that with a few per cent negative correction es to ^^(Mg) we can ac- 
commodate lower values of as{Mz) with no problem. The best fit values of e^ are presented 
in figures |lO|a-b. Clearly, as squarks get lighter the SUSY correction to mb{Mz) (fig-i) has 
to be increasingly reduced by lower values of as{Mz) (fig.|^), which in turn requires a more 
substantial departure from the gauge coupling unification (fig. [TO|) in the respective SUSY 
parameter subspace. 



5 Conclusions 

In summary, we have presented the results of the MSSM analysis focusing on the constraint 
imposed by the measured value of the BR{b — > sj). The analysis was motivated by the 
best of the S'O(IO) SUSY GUT models studied previously. Large tan/3 was assumed as a 
consequence of simple SO(IO) fiavor dynamics. We showed that the effective amplitude for 
the inclusive b —>■ s'-f decay can be of either sign in such a scenario. In the case, when the 
sign is the same as in the SM calculation, the often neglected inter-generational CL-ti squark 
mixing and the NLO correction tend to increase the allowed SUSY parameter space. The 
best fits favor light Higgs spectrum since a substantial charged Higgs contribution to b ^ sj 
helps to counterbalance the chargino contribution. In the complementary case, when the 

^"^The sign of the chargino induced correction to irib is determined by the product At^ (the dominant part 
comes from the same diagram as in fig.||c with sl replaced by &l and no photon leg attached). Our best 
fits always run into the region where At{Mz) turns out negative, in the conventions which maintain positive 
gluino mass parameter. 
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sign of the Cj coefficient is ffipped by a large chargino contribution, both the CL-ti mixing 
and the NLO correction work the other way and tend to reduce the allowed SUSY parameter 
space. The charged and pseudoscalar Higgs masses can be either large or small. In each 
case, the chargino (and neutralino) masses are at the experimental limit in the best fits. 
Because of 5mf^^^ > the best fits result in rather low values of as{Mz), which can be 
obtained with a few per cent negative correction to agiMc)- With these guidelines the large 
tan/5 regime of the MSSM constrained by simple 5*0(10) GUTs remains a viable option for 
physics beyond the Standard Model. 
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Figure 1: Contour lines of constant x^ ^^ the best fits of model 4c for fi{Mz) = (a) 80 GeV, 
(6) 160 GeV and (c) 240 GeV. Solid, dashed and dotted curves correspond to x^= 6, 3 and 
1 per 3 degrees of freedom, respectively. 



17 



(a) 



|l=110GeV 



600 

500 

> 400 

CD 

CD 
^^300 

200 

100 



^ MSSM . ^ SM 

Sign C, = + sign C, 



/ / / / 


( 




/ 
/ 

( 0.4 

\ 
\ 
\ 








\ \ 

\ \ 

> \ 


X'=6X^ 






\0.4 / 













500 



1000 
mo [GeV] 



1500 



2000 



(b) 



1^ 


= 110GeV 




sign C 


MSSM . ^ SM 

V = - Sign C, 


600 
















^y''^ ^^ 


500 








-^^^ 


> 400 

CD 
O 






x^^^^ c^^ 


r^ 


^^300 


/ / 




/ 




- 


200 
H r\c\ 


X^=6.(Lj 


VJ.O 


Vl.O Lo.4 


Lo.3 













500 



1000 
mo [GeV] 



1500 



2000 



Figure 2: x^ contour plots in the best fits of model 4c with the MSSM effective amplitude for 
6 ^ S7 of (a) the same (6) the opposite sign as compared to the SM amplitude. [liMz) = HO 
GeV, which — as described in the text — basically corresponds to the minimum value if /i 
freely varied as an input parameter. As indicated, curves correspond to x^ = 6, 3, 1, 0.4, 0.3 
per Sd.o.f., respectively. 
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Figure 4: Contour lines of constant BR{b -^ 57) x 10^ in the best fits assuming the opposite 
signs of the effective amphtude for this process. For better reference, the x^ contour plots of 
figures 2a and 2b are shown in the background of (a) and (b), respectively, as dotted lines. 
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Figure 5: The same as in fig.§ for the contour plot of r^^-* defined in eq.(|^). The effect of 
the doubled minimum in the lower part of figure (b) corresponds to the possibility of having 
heavy (solid black curves) or light (dashed gray curves) Higgs spectrum for M1/2 < 300GeV. 



21 



(a) 



^=110GeV 



600 



^ MSSM . ^ SM 

Sign C7 = + signC7 




1000 1500 

mo [GeV] 



2000 



(b) 



^=110GeV 



600 



^ MSSM . ^ SM 

Sign C7 = - sign C^ 




1000 1500 

m„ [GeV] 



2000 



Figure 6: The same as in fig.^ for the contour plot of r*^*"-* defined in eq.(|^). 
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Figure 7: The same as in fig.^] for the contour plot of r*^*"^^-* defined in eq.(|Tl|) 
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Figure 8: The same as in fig.^ for the contour plot of 5mf^^^ , the SUSY correction to the 
h quark mass mb{Mz). 
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Figure 9: The same as in fig.§ for the contour plot of as{Mz) 
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Figure 10: The same as in fig^ for the contour plot of 63, the GUT scale threshold correction 

to a,. 
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